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Closed-form  and  finite-element  solutions  are  presented  for  thermal  bend¬ 
ing  and  stretching  of  laminated  composite  plates.  The  material  of  each  layer 
is  assumed  to  be  elastically  and  thermoelastically  orthotropic  and  bimodular, 
l.e.,  having  different  properties  depending  upon  whether  the  fiber-direction 
normal  strain  is  tensile  or  compressive.  The  formulations  are  based  on  the 
thermoelastic  version  of  the  Whitney-Pagano  laminated  plate  theory,  which  in¬ 
cludes  thickness  shear  deformations.  Numerical  results  are  obtained  for  de¬ 
flections  and  neutral -surface  positions  associated  with  normal  strains  in  both 
of  the  in-plane  coordinates.  The  closed-form  and  finite-element  results  are 
found  to  be  in  good  agreement.  ^ 

1  INTRODUCTION 


As  the  use  of  fi_ber-reinforced  composite  materials  becomes  more  wide¬ 
spread,  there  is  ever  increasing  interest  in  predicting  the  thermostructural 
behavior  of  plates  constructed  of  such  materials.  One  of  the  interesting 
characteristics  of  certain  fiber-reinforced  composite  materials  is  that  they 
exhibit  quite  different  elastic  properties  when  loaded  along  the  fiber  direc¬ 
tion  in  tension  as  opposed  to  compression.  This  was  observed  experimentally 
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for  cord-rubber  composites  by  Clark  (lj^and  Patel  et  al.  (2).  he  first 
attempt  to  formulate  a  theory  of  elastic  behavior  of  such  materials  was  made 
by  Ambartsumyan  (3),  and  a  comprehensive  theory  consistent  with  experimental 
results  was  presented  in  (4). 

Although  there  have  been  a  number  of  thermoelastic  analyses  of  bimodulus- 
material  plates  (5-9),  they  have  all  been  limited  to  isotropic  bimodulus 
materials  and  temperature  changes  symmetric  about  the  midplane  of  the  plate. 

The  present  work  is  more  general  than  any  existing  analyses  known  to  the 
authors  in  four  different  respects: 

1.  The  material  of  each  layer  is  elastically  orthotropic  and  bimodular, 

2.  The  material  of  each  layer  is  thermoelastically  orthotropic  and  bi¬ 
modular,  as  shown  to  be  possible  physically  in  a  recent  micromechanics 
analysis  (1_0). 

3.  Both  single-layer  orthotropic  and  cross-ply  laminated  plate  con¬ 
structions  are  considered. 

4.  Transverse  shear  deformation  is  included. 

5.  Temperature  changes  through  the  thickness  as  well  as  in  the  plane 
are  considered. 

2  GOVERNING  EQUATIONS 

The  equations  developed  here  constitute  the  thermoelastic  extension  of  the 
Whitney-Pagano  shear-flexible  laminated  plate  theory  (VI)  . 

The  origin  of  a  Cartesian  coordinate  system  is  taken  to  be  in  the  mid¬ 
plane  (xy  plane)  of  the  plate  with  the  z  axis  being  normal  to  this  plane  and 
directed  downward. 

The  thermoelastic  constitutive  relations  for  each  layer  (z)  are  taken  to 
f  References  are  given  in  Appendix  A. 


be  orthotropic  and  bimodular  as  follows: 
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Here  a  and  a  are  In-plane  normal  stresses,  t  Is  In-plane  shear  stress, 
x  y  xy 

tyZ  and  tx2  are  transverse  shear  stresses,  ex  and  ey  are  In-plane  normal 
strains,  yxy  is  in-plane  engineering  shear  strain,  yyz  and  yxz  are  transverse 
engineering  shear  strains,  T  is  temperature  change  measured  from  the  strain- 
free  temperature,  the  C's  are  Cauchy  elastic  shear  stiffnesses,  the  Q's  are 
plane-stress  reduced  stiffnesses,  and  the  a's  are  thermal -expansion  coeffi¬ 
cients.  The  third  subscript  in  Q-jj^  and  (and  second  in  a^)  refers  to 

the  bimodular  characteristics:  k*l  denotes  properties  associated  with  fiber- 
direction  tension,  k*2  denotes  fiber-direction  compression,  and  k»0  signifies 
that  the  property  does  not  depend  upon  fiber-direction  strain.  Also,  sub¬ 
script  t  refers  to  the  layer  number,  i.e.,  2=1  and  2  for  a  two-layer  laminate. 

The  stress  resultants  and  stress  couples  and  thermal  stress  resultants 
and  couples  are  defined  in  the  standard  manner,  i.e.. 


fh/2 

(VWVV  *  J_h/2  (ax’VTxy,Tyz,Txz)dz 


rh/2 

(VW  *  J_h/2  ^°x,0y’Txy^2 
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t  t  rh/2 

(Nx’Mx>  "  i-h/2  +  ^i2ki*2lcA  ^ *z  ^  dz  (•*) 

T  T  fh/2 

(Ny*My)  -  J_h^2  (Q12kA°lkA  +  ^22k2°2kA^1,Z^T  dz 

where  h  is  the  total  laminate  thickness. 

The  displacements  u,v,w  In  directions  x,y,z  at  an  arbitrary  location 
(x,y,z)  are  given  by 

u(x,y,z)  *  u°(x,y)  +  Zi|<x(x,y) 

v(x,y,z)  *  v°(x,y)  +  z<|>y(x,y)  (4) 

w(x,y,z)  =  w°(x,y) 

Here  u°,v°,w0  are  the  displacements  at  the  midplane,  and  a  and  *  are  the 

x  y 

bending  slopes. 

Substituting  equations  (4)  into  the  well-known  linear  strain-displacement 
relations  of  elasticity  theory  in  Cartesian  coordinates,  and  then  using 
equations  (1)  and  (2)  in  equations  (3),  one  obtains  the  thermoelastic  con¬ 
stitutive  relations  for  the  laminate  as  follows: 


f  T  ^ 

Nx  +  Ni 

An 

A12  0 

Bu/h 

B12/h 

0 

' 

V"J 

Au 

A22  0 

B12/h 

B22/h 

0 

Nxy 

0 

0  a66 

0 

0 

B66/h 

vVu°y 

<  T 

s  = 

< 

V 

(MX  +  Mx)/h 

Bn/h 

Biz/h  0 

Dn/h2 

Di2/h2 

0 

hip 

x,x 

(My  +  Mj)/h 

Bu/h 

B22/h  0 

DI2/h2 

D22/h2 

0 

hip 

y,y 

Vh 

v  xy  > 

0 

0  B66/h 

0 

0 

B66/h2 

hu>  +  hip 

yix  w**y) 

4 


Here  a  comma  denotes  differentiation,  i.e.,  (  )>x  =  a(  )/3x. 

As  usual,  the  stretching,  bending-stretching  coupling,  bending  and  trans¬ 
verse  shear  stiffnesses  for  the  laminate  are  defined  as 

fh/2 

(A10*Bij’Dij»Sii)  *  J_h/2  ^ij’Qijz,Qijz2,KiCii^dz  ^ 

where  Kz  is  the  transverse  shear  coefficient  to  provide  for  the  nonuniform 
transverse-shear  strain  distribution  through  the  thickness  (12,). 

Neglecting  distributed  forces,  body  moments,  and  inertial  effects,  one 
can  write  the  equations  of  equilibrium  as 


N  +  N  *0;N  +  N  *  0  ;  Q  +0  *0 

x,x  xy,y  *  xy,x  y,y  ’  yx,x  vy,y 


M  +  M  -  Q  »  0  ;  M  +  M  -0*0 

x,x  xy,y  yx  *  xy,x  y,y  yy 


Substituting  equations  (5)  and  (6)  into  equations  (8),  one  obtains  the 
following  coupled  thermoelastic  equilibrium  equations  in  terms  of  the 
generalized  displacements: 


,  (k,i*1  ,2, 3, 4, 5) 


(9) 


where  [L^]  Is  a  symmetrlx  linear  differential  operator  matrix  with  the 
following  elements: 


L11  2  AlldX  +  A66dy 

l12  =  (A12  +  A66)  dxdy 
L13  5  0 

L14  =  (B12  +  Bgg)  dxdy 
L15  =  Budx  +  B66dy 

1-22  5  Aggd*  +  A22dy 

L23  5  0 

1-24  -  B$$dx  +  B22dy 


L25  2  1-14 

L33  2  -  S55d^  -  S44d* 

L34  =  -  S44dy 

L35  2  -  Ss5dx 

2  2 

1-44  =  Dggdx  +  022dy  "  $44 
1-45  2  (0l2  +  Ogg  )  dxdy 
L55  =  Dlldx  +  Dggdy  "  $55 
dv  =  3/3X  . 


(10) 


Except  for  the  presence  of  the  generalized  thermal-force  terms  appearing 
on  the  right  side  and  the  absence  of  mechanical  pressure,  equations  (9)  are 
identical  to  those  presented  recently  in  (13 ). 


3  CLOSED-FORM  SOLUTION 

Guided  by  the  closed-form  solution  presented  in  (1^)  for  a  freely 
supported,  laminated,  bimodulus,  rectangular  plate  subjected  to  a  sinusoidally 
distributed  normal  pressure,  we  consider  here  the  same  class  of  plate  and 
geometry  but  subjected  to  sinusoidal  distributions  of  midplane  temperature  and 
temperature  gradient  through  the  thickness.  Thus,  the  temperature  distribution 
Is  given  by 

T(x,y,z)  =  T0(x,y)  +  (z/hjT^x.y)  (11) 

where,  here 
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(12) 


Tq  »  Tq  sin  ax  sin  sy,  Tx  *  sin  ax  sin  Sy 
«■*/»,  6  *  ir/b 

and  a  and  b  are  the  lengths  of  the  plate  sides  parallel  to  the  x  and  y  axes, 
respectively. 

The  boundary  conditions  are  freely  supported,  i.e.,  simply  supported 
flexurally  and  unrestrained  in  the  in-plane  directions  normal  to  the  edges. 
Along  the  edges  located  at  x*0  and  x*a: 

Nx  *  v°  "  w  *  Mx  *  ^y  *  0  O3) 

Along  the  edges  located  at  y*0  and  y*b: 

u°  ■  N  »  w  ■  M  ■  *  0  (14) 

For  the  temperature  distribution  and  boundary  conditions  specified  above, 
the  governing  equations  are  satisfied  exactly  in  closed  form  by  the  following 
distributions  of  generalized  displacements: 

u°  *  U  cos  ox  sin  Sy 

v°  *  V  sin  ox  cos  ey 

w  *  W  sin  ox  sin  sy  (15) 

hi|iy  »  Y  sin  ox  cos  sy 

h<|ix  *  X  cos  ox  cos  sy 

The  values  of  the  coefficients  U,V,W,Y,  and  X  are  obtained  by  solving  the 
following  linear  matrix  equation  obtained  by  substituting  solutions  (15)  into 
governing  equations  (9): 
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where  CC^_ ]  is  a  symmetric  matrix  containing  the  following  elements: 


'k  i 

Cn  -  -  Auo2  -  Agge2 
C12  -  "  (Aj2  +  Agg)aB 
C13  =  0 

C14  =  -  (B12  +  Bgg)oB 
Ci 5  2  ~  Bna2  -  B66e2 

C22  -  “  Aggo2  -  A22B2 

C23  =  0 

C24  2  "  Bgga2  -  B2262 


C25  -  ^14 

C33  2  -  Sgsa2  -  S4482 

C34  2  “  S448 

C35  2  “  S5502 

C44  =  -  Dgga2  -  O2262  “  S44 

C45  =  -  ( D i 2  +  Dgg)o8 

C55  2  "  Dll®2  -  DggB2  -  S5S 


(16) 


07) 


and  the  quantities  on  the  right  side  are  defined  by 


nI  v  *  nJ  /cos  ox  sin  ey  ,  nJ  *  nJ  /sin  ox  cos  ey 

A, A  AjA  Jr  tjf  J  *Jr 

fil  v  =  mJ  ¥/cos  ox  sin  ey  ,  u  *  /sin  ox  cos  ey 

A*A  AjA  J  9J  J  »jr 


(18) 


Even  in  the  case  of  a  single-layer  plate  of  bimodulus  material  as  con¬ 
sidered  in  (14),  the  plate  stiffnesses  A^.,  B^.,  and  D^.  depend  upon  the 
neutral -surface  position  (denoted  by  zn)  associated  with  the  fiber-direction 
strain.  In  the  case  of  a  cross-ply  plate,  as  considered  in  (13),  these 
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stiffnesses  depend  upon  both  znx  (In  the  layer  oriented  in  the  x  direction) 
and  zny  (In  the  other  layer,  which  is  oriented  in  the  y  direction).  It  is 
noted  that  the  transverse  shear  stiffnesses  are  unaffected  by  zpx  and  zfly. 

To  determine  the  neutral -surface  location,  we  use  the  kinematics  of 
deformation,  l.e., 

«X  *  0  ■  U°X  *  2*X,X  '•  «y  *  0  •  »?y  *  **y,y  (’9) 


Then,  using  equations  (15),  we  obtain 


nx 


uV*x,x  ■  *  hU'X 


z  =  -  v,  U  =  -  hV/Y 
ny  ’y/vy,y 


(20) 


In  the  present  thermoelastic  case,  not  only  the  A..,  8..,  and  0-x  stiff- 

I J  I J  I J 

nesses  depend  upon  znx  and  zny,  but  also  do  the  generalized  thermal  forces 
nT  and  m|.  The  detailed  form  of  these  dependences  are  a  function  of  the 
nature  of  the  signs  of  znx  and  zpy;  they  are  developed  for  various  possible 
cases  in  Appendix  B. 

In  principle,  one  could  develop  a  set  of  explicit  relations  for  deter¬ 
mining  zpx  and  zpy.  However,  the  extreme  complexity  of  the  algebraic  struc¬ 
ture  of  the  resulting  equation  renders  this  approach  impractical.  Thus,  an 
Iterative  procedure  analogous  to  that  used  in  (12,13)  for  mechanical  loading 
is  adopted  here.  The  procedure  turns  out  to  be  computationally  quite 
efficient. 


4  FINITE-ELEMENT  FORMULATION 

As  pointed  out  in  the  previous  section,  an  exact  closed-form  solution 
to  equations  (9)  can  be  obtained  only  under  special  conditions  of  geometry, 
edge  conditions,  loadings,  and  lamination.  Here  we  present  a  simple  finite- 
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element  formulation  which  does  not  have  any  limitations  (except  for  those 
Implied  In  the  formulation  of  the  governing  equations) (J5.,J6). 

Suppose  that  the  region  "R  Is  subdivided  Into  a  finite  number  N 

of  subregions:  finite  elements,  "Re  (e=l ,2, . . . ,N) .  Over  each  element  the 

generalized  displacements  (u0,v0,w,4ix,i py)  are  interpolated  according  to 

r  i  r  j  s 

u°  *  Z  ,  v°  *  I  v.$.  ,  w°=  z  w.<|^ 

i  i  1 

P  3  P  3 

K  =  5  =  2  ey.*. 


(21) 


where  4>“  (o*l,2,3)  is  the  interpolation  function  corresponding  to  the  i-th 
node  in  the  element.  Note  that  the  in-plane  displacements,  the  transverse 
displacement,  and  the  slope  functions  are  approximated  by  different  sets  of 
interpolation  functions.  While  this  generality  is  included  in  the  formu¬ 
lation  (to  indicate  the  fact  that  such  independent  approximations  are  possible), 
we  dispense  with  it  in  the  interest  of  simplicity  when  the  element  is  actually 
progranmed  and  take  <f>.j  8  $.j  «  (ras=p).  Here  r,s,  and  p  denote  the  number 
of  degrees  of  freedom  for  each  variable.  That  is,  the  total  number  of  degrees 
of  freedom  per  element  is  2r  +  s  +  2p. 

Substituting  equation  (21)  into  the  Galerkin  integrals  associated  with 
the  operator  equation  (9),  which  must  also  hold  in  each  element  He, 

([L]{6>  -  {f }){$}dxdy  =  0  (22) 

and  using  integration  by  parts  once  (to  distribute  the  differentiation  equally 
between  the  terms  in  each  expression),  we  obtain 
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[K11][K12][K13][K1‘t][K15]”|  {u°>|  T  {F1}' 

[K22][IC23][K2^][K25]  {v°}  {F2} 

[K33][K34][K35]  ,  {w°}l  *  \  {F3}  r  ('< 

Symmetric  [K44][K45]  {F4} 

CK”]Je  [(FS>Je 

where  the  {u0},  {v0},  etc.  denote  the  columns  of  the  nodal  values  of  u°,  v°, 

respectively,  and  the  elements  K?!  (a,esl,2 . 5)  of  the  symmetric  stiffness 

matrix  and  F®  of  the  force  vector  are  given  by 


*13  ■  flnGU  *  Mj 
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Kl?  *  0 
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where 


G?3  ■  j„  dx  dy 
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(26) 


( i-1 ,2, . . . >s  »  j*l»29...jt) 


XX  X  1  2  3 

and  G^ .  «  G^,  etc.  In  the  special  case  in  which  <j>^  =  <j>.  *  all  of  the 
matrices  in  equation  (26)  coincide. 

In  the  present  study,  elements  of  the  serendipity  family  are  employed 
with  the  same  interpolation  for  all  of  the  variables.  The  resulting  stiffness 
matrices  are  20  by  20  for  this  four-node  element  and  40  by  40  for  the  eight- 
node  element.  Reduced  integration  (V7,l_8)  must  be  used  to  evaluate  the  matrix 
coefficients  in  equation  (24).  That  is,  if  the  four-node  rectangular  element 
is  used,  the  lxl  Gauss  rule  must  be  used  in  place  of  the  standard  2x2  Gauss 
rule  to  numerically  evaluate  the  coefficients  K... 

*  J 

Substituting  solution  (23)  into  equations  (20),  we  get 


nx 
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5  NUMERICAL  RESULTS 

As  a  check  on  the  validity  of  the  equations  and  their  computational 
implementation,  it  is  desirable  to  compare  the  present  predictions  with  those 
given  in  the  literature  for  special  cases.  Apparently  there  is  a  dearth  of 
solutions  of  plates  bent  by  a  sinusoidally  distributed  thermal  gradient.  How¬ 
ever,  it  was  possible  to  compare  with  isotropic,  thin  plate  results  given  by 
Boley  and  Weiner  (19.)  as  listed  in  Table  1.  As  can  be  seen,  the  agreement  is 
quite  good. 

As  examples  of  orthotropic  bimodulus  materials,  the  same  materials  as 
considered  in  (13, 1£)  are  used, namely,  aramid-rubber  and  polyester-rubber. 

The  elastic  properties,  obtained  from  experimental  results  of  (2.)  using  the 
data-reduction  procedure  presented  in  (4),  are  listed  in  Table  2.  Unfortunately 
there  are  no  measured  values  available  for  the  thermal -expansion  coefficients 
of  these  materials.  However,  the  micromechanics  analysis  of  bimodular  action 
presented  in  (10.)  suggests  that  the  thermal -expans ion  coefficients  of  these 
materials  should  also  depend  upon  the  sign  of  the  fiber-direction  strain.  Thus, 
for  the  numerical  calculations  presented  here,  it  is  assumed  that  the  following 
relationships  hold: 

a^/oj  *  0.5  ;  z\/a\  *  1.0  ;  a\/a\  *  0.1 

Numerical  results  are  presented  for  relatively  thick  plates  (b/h*10) 
with  a  temperature  distribution  having  a  temperature  gradient  through  the 
thickness  but  no  mean  temperature  change,  i.e.,  fQ*0.  Table  3  gives  results 
for  single- layer  orthotropic  plates  of  various  aspect  ratios. 

As  can  be  seen,  there  is  good  agreement  between  both  dimensionless  de¬ 
flections  and  neutral -surface  locations  as  predicted  by  the  finite-element 
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and  closed-form  solutions.  In  fact,  the  largest  difference  appearing  in  the 
table  are  0.1 5%  for  W,  0.18  for  Z  ,  and  1.3*  for  Z  . 

*  y 

Table  4  shows  the  finite-element  and  closed-form  solutions  for  dimension¬ 
less  deflections  and  neutral -surface  locations  of  two-layer  cross-ply  (0°/90°) 
simply-supported  rectangular  plates  subjected  to  sinusoidal  thermal  loading. 
The  finite-element  results  are  in  close  agreement  with  the  closed-form  results. 

It  should  be  mentioned  that  composite  materials  typically  have  much  lower 
ratios  of  thickness  shear  moduli  to  in-plane  Young's  moduli  than  homogeneous, 
isotropic  materials  (G/E  =  1/3  to  1/2).  In  contrast,  aramid-rubber  has 
Gxz/Ex  *  0.001  when  the  fiber-direction  strain  is  tensile,  but  0.416  for  com¬ 
pressive  fiber-direction  strain. 

Figure  1  presents  the  influence  of  the  aspect  ratio  and  side-to-thickness 
ratio  for  single-layer  and  two-layer  cross-ply  plates  under  sinusoidal  thermal 
loading  (material:  polyester-rubber) .  The  effect  of  thickness  on  the  deflec¬ 
tion  is  more  pronounced  than  the  effect  of  the  aspect  ratio. 

Figures  2  and  3  show  the  effect  of  the  aspect  ratio  (a/b)  and  side-to- 
thickness  ratio  (b/t),  respectively,  on  the  location  of  neutral  surfaces  for 
a  single-layer,  orthotropic,  bimodulus,  simply  supported  rectangular  plates 
under  sinusoidal  thermal  loading  (material:  aramid-rubber,  To=0.0,  Tx=l .0) - 

Similar  results  are  presented  in  Figures  4  and  5  for  a  two-layer  cross- 
ply  (0°/90°),  rectangular  plate  of  polyester-rubber  under  sinusoidal  thermal 
loading.  Note  from  Figure  5  that  the  neutral -surface  locations  are  virtually 
unchanged  for  side-to-thickness  ratio  greater  than  25. 


6  CONCLUSIONS 


On  the  basis  of  the  excellent  comparisons  with  existing  results  for  homo 
geneous-material  plates  and  the  comparisons  with  a  closed-form  solution  for 
laminated  bimodulus  composite-material  plates,  the  mixed  finite  element  for 
thick  plates  as  presented  here  is  considered  to  be  validated. 

7  ACKNOWLEDGMENTS 

The  authors  are  grateful  to  the  Office  of  Naval  Research,  Structural 
Mechanics  Program  for  financial  support  thrpugh  Contract  N00014-78-C-0647  and 
to  the  University's  Merrick  Computing  Center  for  providing  computing  time. 


APPENDIX  A 
REFERENCES 

(1)  CLARK,  S.K.,  'The  plane  elastic  characteristics  of  cord-rubber  laminates'. 
Textile  Research  J.  1963  33,  295-313. 

(2)  PATEL,  H.P.,  TURNER,  J.L.,  and  WALTER,  J.D.,  'Radial  tire  cord-rubber 
composites'.  Rubber  Chem.  and  Tech.  1976  49,  1095-1110. 

(3)  AMBARTSUMYAN,  S.A.,  'The  basic  equations  and  relations  of  the  different- 
modulus  theory  of  elasticity  of  an  anisotropic  body'.  Mechanics  of  Solids 
1969  4  (No.  3),  48-56. 

(4)  BERT,  C.W.,  'Models  for  fibrous  composites  with  different  properties  in 
tension  and  compression,  J.  Eng.  Matls.  and  Tech.,  Trans.  ASME  1977  99H, 
344-349. 

(5)  AMBARTSUMYAN,  S.A.,  'The  equations  of  temperature  stresses  of  different- 
modulus  materials'.  Mechanics  of  Solids  1968  3,  (No.  5),  58-69. 

(6)  AM8ARTSUMYAN ,  S.A.,  'Equations  of  theory  of  thermal  stresses  in  double¬ 
modulus  materials',  in  BOLEY,  B.A.  (ed.),  Themoinelasticity  (Proc.,  IUTAM 
Sympos.,  E.  Kilbride,  Scotland,  June  25-28,  1968)  1970,  17-32,  Springer- 
Verlag,  Wien. 

(7)  KAMIYA,  N.,  'Thermal  stress  in  a  bimodulus  thin  plate'.  Bulletin  de  V 

academic  Polonaise  des  Sciences ,  Serie  des  sciences  techniques  1976  24, 
365-372.  ~ 

(8)  KAMIYA,  N.,  'Energy  formulae  of  bimodulus  material  in  thermal  field',. 
Fibre  Sai.  and  Tech.  1978  11,  229-235. 


15 


1 


(9)  KAMIYA,  N.,  'Bimodulus  thermoelasticity  considering  temperature-dependent 
material  properties',  BERT,  C.W.  (ed.)»  Mechanics  of  Bimodulus  Materials 
Dec.  2-7,  1979,  29-37,  ASME  Winter  Annual  Meeting,  New  York. 

(10)  BERT,  C.W.,  'Micromechahics  of  the  different  elastic  behavior  of  filamentary 
composite  in  tension  and  compression',  in  BERT,  C.W.  (ed.),  Mechanics  of 
Bimodulu8  Materials  Dec.  2-7,  1979,  17-28,  ASME  Winter  Annual  Meeting, 

New  York. 

(11)  WHITNEY,  J.M.  and  PAGANO,  N.J.,  'Shear  deformation  in  heterogeneous  aniso¬ 
tropic  plates',  J.  Appl.  Meah.  1970  37_,  1031-1036. 

(12)  WHITNEY,  J.M. ,  'Shear  correction  factors  for  orthotropic  laminates  under 
static  loads',  J.  Appl.  Meah.  1973  40,  302-304. 

(13)  BERT,  C.W.,  REDDY,  J.N.,  REDDY,  V.S.,  and  CHAO,  W.C.,  'Analysis  of  thick 
rectangular  plates  laminated  of  bimodulus  composite  materials',  AIAA/ASME/ 
ASCE/AHS  2lst  Structures,  Structural  Dynamics  and  Materials  Conference, 
Seattle,  May  12-14,  1980. 

(14)  BERT,  C.W.  and  KINCANN0N,  S.K.,  'Bending-extensional  coupling  in  elliptic 
plates  of  orthotropic  bimodulus  material'.  Developments  in  Mechanics,  U) 
(Proc.,  16th  Midwestern  Mechanics  Conference,  Sept.  19-21,  1979),  7-11, 
Kansas  State  Univ,,  Manhattan,  Kansas. 

(15)  RED0Y,  J.N.,  'A  penalty  plate-bending  element  for  the  analysis  of  laminated 
anisotropic  composite  plates',  Int.  J.  Burner.  Meth.  Engng.  1980,  in  press. 

(16)  REDDY,  O.N.  and  HSU,  Y.S.,  'Effects  of  shear  deformation  and  anisotropy 
on  the  thermal  response  of  layered  composite  plates',  j.  Thermal  Stresses , 
to  appear. 

(17)  ZIENKIEWICZ,  O.C.,  TAYLOR,  R.L.,  and  TOO,  J.M.,  'Reduced  integration 
technique  in  general  analysis  of  plates  and  shells',  Int.  J.  Burner.  Meth. 
Engng.  1971  3,  575-586. 

(18)  REDDY,  J.N.,  'A  comparison  of  closed-form  and  finite-element  solutions  of 
thick,  laminated,  anisotropic  rectangular  plates'.  Report  OU-AMNE-79-19, 
University  of  Oklahoma,  Norman,  OK,  Dec.  1979. 

(19)  B0LEY,  B.A.  and  WEINER,  J.H.,  Theory  of  Thermal  Stresses ,  Wiley,  New  York, 
1960,  389-391. 


APPENDIX  B 

DERIVATION  OF  EXPRESSIONS  FOR  THERMAL  FORCES  AND  MOMENTS 


Case  I 


For  Case  I,  zpx  >  0  and  zn^  <  0  with  znx  governing  layer  1  (0°)  and  zny 
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layer  2  (90°), 
z. 


T  f  n*  , 

"x  *  (Q 
x  J-h/2 


rZnx 


'1122  b122  +  Q1222  °222)T  ***]  (Q1112  a112  +  Q1212  a212)T  dz 


r° 

,)T  +  I 


ny 

h/2 

znx 


(8-1 ) 


ln  rh/2 

(Q1 121  *121+Q1221  °221  T  dz  +  J  (Q1111  °m+Q1211  a211)T  dz 


Let 


(QH22  °122  +  Q1222  a222)  *  6122  *  (Q1112  °112  +  Q1212  a212)  "  B112 

(B-2 

(Q1121  °121+Q1221  a221}  "  S1 21  *  (Q1111  0,111+Q1211  a211}  "  8111  •  etc- 
Then, 

Nx  *  [“122  To<2ny  *  */2>  *  *112  V°-2ny>  *  *121  V2„x-°> 

*  *111  To(h/2-2nx>  *  s122(Tl/2H»2ny- h2/4> 

*  *„2(V2h»°-2^)  *  »121<V2hH4-0) 

+  s'|i-|(Ti/2h) (h2/4 -  zzx)]sin  ax  sin  By 
or 

Nx  *  C(*122**lll>(Toh/2)  *  <e121  *  Blll>  Vnx  *  (*122-*112)  Vny 

*  <B111  *  S,22)<T,h/8>  *  <*,2, -*,„)<  VnV2h>  (B-3) 

+  ^B122  ”  ®H2^^1zny^^^s^n  oX  Sy 
Similarly, 
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Ny  "  ^B222  +  B211^Toh/2)  +  ^22T6211)  Vnx  +  (6222’6212)  Vny 


+  (®211  "  ®222^Tlh/8^  +  ^B221  "  B211^TlZnx/2h^  +  ^B222"6212^ 


(Tlzny/2^)^sin  aX  sin  By 


T  fh/2  f0  r  nx  h/2 

\  "  J _  B122  Tz  dz  +  Jz  B112Tzdz+J0  6121Tzdz+  z  Blll  Tz  dz 


^B111  '  8i22)(Toh2/8)  +  (b121  ■  Bm)(Toznx/2)  +  (b122  '  B112)(ToZny/2) 


+  (8122  +  8m)(Tih  /24)  +  (b121'  6111)(Tlznx/3h) 


+  (8i22"  B112^Tlzny/3h^sin  oX  sin  By 


Similarly, 


My  *  ^®211  "  B222^Toh2/8)  +  *B221  "  B211  ^ToZnx/2*  +  ^s222  "  B2l2^ToZny/2^ 


+  ^222  *  ®2U^1^  +  (8991  *  89yj)(T«|  zrtv/3h) 


221  B2]]M  1  nx' 


^  ( 8222  ”  ®2] 2} (^1 z ny^3h) 3sin  ax  cos  By 


Using  the  above  equations  In  conjunction  with  equations  (12)  and  (18),  we 


obtain  the  following: 


C  “  °^6122  +  B111)(Toh/2)  +  (B121  "  B111)Toznx  +  (s122  “  B112)Tozny 

+  (B111  “  B122)(Tlh/8)  +  (B121  '  Blll^Vnx/2h) 

+  C6122  “  B112)(Tlzny/2h^ 


I 


(B-4) 


(B-5) 


(B-6) 


(B-7) 


+  >1  *  &K&222  +  0211^o^2^  +  ^8"i  “  )^„z„v  +  (B999  "  891 9)Tftz 


221  211 '  0  nx  '222  B212'Vny 


+  (8211  -  6222)(T1h/8)  +  (622i  *  s211^Vnx/2h* 

+  (8222-S212)(TizJy/2h)] 

«I  *  *  aC(0-|ii  “  8-j 22) (T0h2/8)  +  (8191  -  8ni)(TAZ^v/2)  +  (8i99-8n9) 


(B-8) 


121  111M  0  nx' 


122  °112' 


*Vny/2*  +  ^122  +  8111^TlhZ/24^  +  ^s121  "  8111  ^Vnx/3h) 

+  (8l22'6112)(Tlzny/3h^ 

*y,y  *  8^8211  "  6222^V2/8)  +  (s221  “  s211  ,(Toznx/2^  +  ^6222  '  s212^ Vny/2) 


( B— 9 ) 


(B-10) 


+  (»222  ■*"  S21 1  )(Tlh  /24)  +  (S221  -  S2iiMVnx/3h)  +  ^s222  “  s212^Tlzny/3h^ 

In  a  similar  way,  one  can  obtain  the  expressions  for  the  above-mentioned 
quantities  for  the  remaining  seven  cases  as  follows: 

fiaaJi  unx>0>  V0) 


C  *  at(s122  +  0lll)(Toh/2)  +  (0121  *  8lll)(Vnx) 

+  ^8111  ‘  fl122^Tlh/8)  +  (s121  ‘  slll^Vnx/2h)3 

Ry,y  *  C ^222  +  a211)(Toh/2)  *  (S22i  *82u)(Vnx) 

+  (s211  *  8222)(V/8)  +  (0221  “  8211^Vnx/2h^ 

lx  *  ®^B111  ”  8122)(T0h2/8)  +  (B-,21  ’  8mHToznx/2) 

+  (8122  +  8in)(T1h2/24)  ♦  (8i2T  slll)(Tlznx/3h^ 

fty,y  *  8^8211  “  8222^ty2/8)  +  (0221  “  8211^Vnx/2^ 

+  (S222- S211)(T1h2/24)  +  (S221  -  ®2n)^izJx/3h)3 


(8-11) 
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Case  ^  (znx<0’  zny>0) 

Rx,x  *  °^8122  +  8111)(Toh/2)  +  ^a121  "  8lll^Vny* 

+  (b122 "  8112^Vnx^  +  (Sin  ’  »i22^Tlh/8^ 

+  (B12i  "  Sni)(Tizny/2h)  +  (®122"  8112^T1znx/2h^ 

N^y  *  sC( 6222  +  6211  HTQh/2)  +  (s221  -  B211)(t0*ny* 

+  (S222  "  a212^Vnx^  +  (8211  "  8222^V/8* 

+  (0221  "  621l)(Vny/2h)  +  (s222  *  8212^1znx/2h^ 
fix,x  "  “C(0in  "  0i22^^oh2/8^  +  ^e121  "  alll^Vny/2* 

+  (0122  “  s112^Vnx/2^  +  (®122  +  Slll^V^24* 

+  (0]21  ’  a171^TlZny/3h^  +  ( s122  “  81 12^1Znx/3h^ 

”y,y  *  ^S211  "  8222^V*2/8^  +  ^  s221  "  6211^ToZny/2* 

+  (0222"  a212^ToZnx/2^  +  ^ e222  +  8211^Tlh  /24^ 

+  (®221  “  8211^TlZny/3h)  +  (8222 ' 8212^1Znx/3h^ 

Case  IV  (znx<0,  zny<0) 

C  *  a^8122  +  8lll^Toh/2*  +  *8122"  8121^Vny* 

+  (0111  "  0l22)(Tlh/8)  +  (8122  ”  8121  ^1zny^2^ 

Ry,y  "  b^8222  +  8211)(Toh/2)  +  (®222"  8221^Vny* 

*  ( 821 1  *  8222 ^ (^l ^8)  +  ( 8222  ”  8221^1zny^2^ 


(B-12) 
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x  *  °[(0m  "  8122^oh2^  +  ^8122’  8121^ozny^ 

+  ^8122  +  flm,(Tlh2/24)  *  ^8122  “  8121  (B“13) 

^y»y  *  ®^( S21 1  "  8222^oh2>^  +  ^8222*  8221^oZny^ 

+  (®222  +  8211^1^2^24^  +  ^8222  ”  8221  ^1zny^8^ 

For  neutral  surface  going  out  of  plane. 

Case  V  (znx>0.5,  zny<"0,5^ 

Rx,x  *  °C(0121  +  8112^V2*  +  ^8121  *  8112^V8^ 

*J,y  ’  8^8221  +  8212)(V2)  +  (8221  ’  8212)(V8)]  (B.14) 

fix,x  *  °^8121  '  S112^V8^  +  ^8121  +  8112^V24^ 

fly  y  *  s£(  ^221  ”  s212)(V8)  +  (8221  +  8212)(T1/24^ 

Case  VI  (2nx<-°-5*  2ny>0-5) 

Nx  x  *  aCCs^i  *  fl112)(V2*  +  {S121  '8112)(V8)J 

^y,y  *  8^8221  +  8212^V2^  +  ^8221  *  S212^V8^  (B>15; 

flx, x  *  e‘^8121  “  8112^V8^  +  ^121  +  8112^V24^ 

fly, y  *  8^8221  *  8212^V8^  +  ^8221  +  ®212^V24^ 

Case  VII  (2nx>0.5,  2^0*5) 

flx>x  *  <»C( S-| -j i  +  8H2^^o^2^  +  ^8111  ’  8112^^1^8^ 

Sy,y  *  sC(8211  +  8212)(To/2)  +  (8211  ’8212^V8)]  (B_16 

*  ®C(8ni  ■  ®112^V8^  +  ^81U  +  8112^Tl/24^ 

^y.y  *  8^8211  '  8212^V8^  +  ^82n  +  S212^V24^ 
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Case  VHI  (znx<-0.5,  zny<-0.5) 

*x,x  =  °C(S121  +  B122^V2)  +  (6121  "  b122)(T1/8^ 

*  b[(b221  *  B222)(V2>  *  (b221-S222HT,/8)] 

fix,x  ‘  “C(6121‘8122)<y8)  *  (S,j,  ♦  8I22)(T,/24)J 

Vy  '  s^(6221  '  B222,(V8)  *  (b221  +  a2Z2* <T1/Z4)3 

For  a  single  layer,  change  S))2  to  8,„,  S)22  to  »gt  to  a2„ 

S222  t0  S22T 
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Table  1  Comparison  with  Boley  and  Weiner's  Results  [17]  for  an  Isotropic 
Single-Layer  Thin  Rectangular  Plate  at  Different  Aspect  Ratios 
(E11/E22  *  1-00  »  vi2  *  v2i  *  0*3  .  b/h  »  10) 


Aspect 

Ratio, 

a/b 

Deflection,  W/h  (m! 

*n»1 ) 

Boley  and  Weiner 

Present 

0.5 

0.5300 

0.5264 

1.0 

6.5858 

6.5789 

1.5 

6.3112 

6.3063 

2.0 

2.1104 

2.1058 

Table  2  Elastic  Properties  for  Two  Tire-Cord/Rubber,  Unidirectional ,  Bimodulus 
Composite  Materials® 


Aramid- 

■Rubber 

Polyester-Rubber 

Property  and  Units 

k=l 

k-2 

k*l 

k*2 

Longitudinal  Young's  modulus,  GPa 

3.58 

0.0120 

0.617 

0.0369 

Transverse  Young's  modulus,  GPa 

0.00909 

0.0120 

0.00800 

0.0106 

Major  Poisson's  ratio,  dimensionless^ 

0.416 

0.205 

0.475 

0.185 

Longitudinal -transverse  shear  modulus,  GPa® 

0.00370 

0.00370 

0.00262 

0.00267 

Transverse-thickness  shear  modulus,  GPa 

0.00290 

0.00499 

0.00233 

0.00475 

aFiber-d1 recti on  tension  is  denoted  by  k*l,  and  fiber-direction  compression 
sby  k-2. 

It  Is  assumed  that  the  minor  Poisson's  ratio  is  given  by  the  reciprocal 
relation. 

°It  Is  assumed  that  the  longitudinal-thickness  shear  modulus  is  equal  to  this 
one. 
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Table  3  Neutral -Surface  Positions  and  Dimensionless  Deflections  for  Rectangular 
Plates  of  Single-Layer  0°  Aramid-Rubber  and  Polyester  Rubber  Determined 
by  Two  Different  Methods  (b/h=10.0,  T^l.O,  To=0.0) 


Aspect 

Ratio 

a/b 

Q* 

Zx 

h 

C.F.** 

** 

F.E. 

★  ★ 

C.F. 

★* 

F.E. 

★★ 

C.F. 

★  ★ 

F.E. 

Aramid 

-Rubber: 

0.5 

0.03478 

0.03486 

0.03059 

0.03059 

0.2853 

0.2807 

0.75 

0.08147 

0.08159 

0.04449 

0.04441 

0.2143 

0.2120 

1.0 

0.1522 

0.1523 

0.06259 

0.06259 

0.1332 

0.1323 

1.25 

0.2485 

0.2485 

0.07970 

0.07968 

0.07699 

0.07699 

1.50 

0.3624 

0.3624 

0.08850 

0.08849 

0.04492 

0.04482 

1.75 

0.4792 

0.4790 

0.08612 

0.08612 

0.02843 

0.02838 

2.0 

0.5880 

0.5876 

0.07491 

0.07491 

0.01987 

0.01986 

Polyester-Rubber: 

0.5 

0.04815 

0.04823 

0.1031 

0.1030 

0.1851 

0.1827 

0.75 

0.1157 

0.1158 

0.1184 

0.1183 

0.1068 

0.1062 

1.0 

0.2160 

0.2161 

0.1308 

0.1308 

0.05299 

0.05276 

1.25 

0.3410 

0.3409 

0.1360 

0.1360 

0.02552 

0.02546 

1.5 

0.4737 

0.4736 

0.1332 

0.1331 

0.01285 

0.01282 

1.75 

0.5975 

0.5970 

0.1234 

0.1233 

0.007329 

0.007329 

2.0 

0.7024 

0.7017 

0.1078 

0.1078 

0.005110 

0.005109 

**W  =  Wh/aftib2 

C.F.  denotes  closed-form  solution;  F.E.  signifies  finite-element  solution 


Table  4  Neutral -Surface  Positions  and  Dimensionless  Deflections  for  Rectangular 
Plates  of  Two-Layer  (0°/90°)  Aramid-Rubber  and  Polyester-Rubber  Under 
Sinusoidal  Thermal  Loading,  as  Determined  by  Two  Different  Methods 
(b/h  »  10.0,  T i  *  1.0,  T0  =  0.0) 


*5)  *  Wh/a^b2 

C.F.  denotes  closed-form  solution;  F.E.  signifies  finite-element  solution 
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cig.  1  Transverse  deflection  vs.  aspect  ratio  (a/b)  and  side-to- 
thickness  ratio  (b/h)  for  single-layer  and  two-layer  cross 
ply  plates  under  sinusoidal  thermal  loading,  (material: 
polyester-rubber) 


Fig.  2  Neutral -surface  location  vs.  aspect  ratio  for  single¬ 
layer  (0°)  rectangular  plates  under  sinusoidal  thermal 
loading,  (material:  aramid-rubber ,  b/h  *  10) 


Fig.  3  Neutral -surface  location  vs.  side-to-thickness  ratio  for 

single-layer  square  plates  under  sinusoidal  thermal  loading 
(material:  aramid-rubber) 
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Neutral -surface  location  vs.  aspect  ratio  for  two-layer 
cross-ply  (0°/90°)  rectangular  plates  under  sinusoidal 
thermal  loading  (material:  polyester-rubber,  b/h  =  10) 


o  □  FEM 
-  CFS 


Neutral -surface  location  vs.  side-to-thickness  ratio  for 
two-layer  cross-ply  {0°/90°)  square  plates  under  sinusoidal 
thermal  loading  (material :  polyester-rubber) 
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finite-element  results  are  found  to  be  In  good  agreement. 
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